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ABSTRACT 


This  report  summarizes  the  major  technical  accomplishments  during 
the  period  from  22  September  1967  to  21  September  1968  under  Contract 
No.  F30602-68-C-0067  (AJRPA  Order  No.  1010).  The  accomplishments  can  be 
grouped  under  three  headings;  namely,  (1)  development  of  a  new  method 
for  directivity  or  signal- to ''noise  optimization  by  spacing  perturbation, 
(2)  development  of  a  new  technique  for  the  optimization  of  arrays  with 
a  large  number  of  elements,  and  (3)  development  of  a  new  method  for 
pattern  synthesis  of  circular  arrays  with  directive  elements.  Essential 
formulations  as  well  as  analytical  justifications  are  outlined. 
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I.  INTRODUCTION 


"Array  Optimization  Criteria"  was  the  title  of  Contract  No. 
F30602-68-C-0067  sponsored  by  the  Advanced  Research  Projects  Agency  under 
ARPA  Order  No.  1010  and  monitored  by  the  Rome  Air  Development  Center  of  the 
Air  Force  System  Command.  This  is  a  technical  report  summarizing  the 
accomplishments  under  the  contract  for  the  period  from  22  September  1967  to 
21  September  1968. 

The  accomplishments  can  be  grouped  under  three  major  tasks.  They  are: 

(1)  Development  of  a  new  method  for  directivity  or  signal-to-noise 
optimization  by  spacing  perturbation, 

(2)  Development  of  a  new  technique  for  the  optimization  of  large  arrays. 

(3)  Development  of  a  new  method  for  pattern  synthesis  of  circular  arrays 
with  directive  elements. 

The  analytical  techniques  developed  for  the  above  tasl  s  are  outlined  in  the 
following  sections. 
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II.  SPACING  PERTURBATION  TECHNIQUES  FOR  ARRAY  OPTIMIZATION1 

Until  the  development  of  the  technique  to  be  outlined  in  this  section, 
array  optimization  usually  starts  with  an  array  of  a  given  configuration, 
and  maximization  of  a  performance  index  is  achieved  by  properly  adjusting 
the  excitation  amplitudes  and  phases  in  the  array  elements.  However,  for 
a  given  set  of  amplitude  and  phase  values,  uniform  spacing  does  not  yield 
the  highest  obtainable  directivity  or  signal-to-noise  ratio.  The  optimum 
element  spacings  can  be  found  by  a  spacing  perturbation  technique.  The 
basis  of  the  technique  lies  in  an  optimization  theorem  which  will  be  proved 
in  the  Appendix.  The  starting  point  can  be  an  array  of  arbitrary  (uniform 
or  nonuniform)  spacings  with  nonoptimum  excitation  amplitudes  and  phases. 
After  an  optimum  set  of  spacings  is  obtained,  the  excitation  amplitudes 
and  phase  shifts  can  be  adjusted  for  further  improvement  in  the  desired 
performance  index;  then  the  optimum  spacings  can  be  recalculated  for  the 
new  excitation  and  the  cycle  repeated  if  desired.  In  particular,  the 
technique  provides  a  method  for  improving  the  directivity  or  the  output 
signal-to-noise  ratio  of  an  array  with  any  given  excitation  and  spacings 
by  spacing  perturbation  until  a  maximum  is  obtained. 

Consider  a  linear  array  of  2N  +  n  identical  antenna  elements  sym¬ 
metrically  located  about  the  origin,  with  n  *  1  when  the  total  number  of 
elements  is  odd,  and  n  =  0  when  the  total  number  of  elements  is  even  and 
the  center  element  is  absent.  Let  6q  be  the  angle  which  the  direction  of 
the  signal  makes  with  the  array  axis.  The  excitation  in  the  mth  element 

from  the  origin  is  I  exp(j$>  ),  where 

m  m 


2 


(1) 


$  =  -  (2-Trd  /X)  cos  0  -  <p 

m  m  o  ib 


In  (1),  is  the  distance  of  the  mth  element  from  the  origin,  X  ia  the 
signal  wavelength,  and  <j>m  is  the  phase  shift  from  th^  cophasal  operation. 
The  array  factor  will  then  be 


where 


N 


E(i|0  =  T~  +  )  I  c°s  (D  ^  "  <f>  )  , 
l  m  mm 

m=l 


ip  =  (2-rrd/X)  (cos  0  -  cos  0q) 


(2) 


(3) 


D  -  d  /d 
m  m 


(4) 


and  d,  a  normalizing  distance,  may  be  any  choice  of  convenience.  For 
example,  if  one  starts  with  a  uniformly  spaced  array,  it  would  be  natural 
to  make  d  the  spacing  between  neighboring  elements.  We  define  the  output 
signal-to-noise  ratio  as  the  ratio  of  the  power  received  per  unit  solid 
angle  in  the  direction  of  the  signal  to  the  average  noise  power  received 
per  unit  solid  angle.  Thus, 


where 
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2  if 


E(0)  | 


TT 


w(0,  <J>) 
g(0,  <p) 
T(0 ,  <j>) 


jjf  d <pf  | E(f|f)  | 2  w(0,  <f>)  sin  0 
0  "0 

-  g(0,  <j>)  T(0 ,  <j>) 

=  element  power  pattern  function 

=  spatial  distribution  function  of 
noise  power 


d0 


(6) 

(7) 

(8) 
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It  is  convenient  to  normalize  g(0,  0)  in  (7)  with  respect  to  its  value 
in  the  direction  of  the  signal;  i.e.,  g(0Q,  $  )  =  1«  The  composite 
function  w(0,  0)  in  (6)  can  be  viewed  as  a  weighting  function  on  the 
array  power  pattern  |  E (ip )  |  2  -  We  note  that  (5)  becomes  the  expression 
for  directivity,  Gq,  when  T(0,  <j>)  =  1;  hence  it  serves  as  the  starting 
point  for  the  optimization  of  both  Gg^  and  Gq.  Clearly,  both  are 
affected  by  the  normalized  element  positions  {D^}  and  the  excitations 


{ 1  »  $  }. 
m  m 


Let  the  perturbed  normalized  element  positions  be  {D^}. 


D  =  D  +  x  , 
m  m  m 


w 


where  Dm  and  x^  represent,  respectively,  the  unperturbed  position  and 
the  spacing  perturbation  for  the  mth  element  and  <<  1.  Substitution 
of  (9)  in  (2)  yields  approximately 


N 

EW  =  E °(«  -  Y  xlIJ  sin  (D%  -  *  )]  ,  (10) 

/  .  mm  mm 

m=l 

where  E°(^)  is  the  unperturbed  array  factor  with  D°  substituted  for 
in  (2).  Using  (10),  we  can  write  (5)  in  the  following  form: 


gsnr 


_ E°(0)!2 

A  -  2x’  F  +  x'  Cx 


(11) 


where 


2tj  it 

I d*  / 

0  0 


E°0(0|2  w(0,  cp)  sin  0  d0 


(12) 
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1»  x2 
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(13) 
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is  the  transpose  of  the  column  matrix  of  spacing  perturbations  x;  3  is  a 
column  matrix  of  typical  element 


m 


/IT  TT 


Im^E°(i|0  w(0,<|>)  sin  (D°tjj 


<j>  )  sin  0  d0  ; 
Tm  * 


(14) 


and  C  °  [c^]  is  an  N  x  N  square  matrix  with 


2tt  4t 


'mk 


h  (  “♦jVk 


^2w(0,<j))  sin  (D^  “  '{'^sinCD^  -  <J>^)sin  0  d0  (15) 


It  can  readily  be  shown  that  C  is  symmetric  and  positive  definite. 
Use  can  then  be  made  of  the  theorem  proved  in  the  Appendix,  which  enables 
us  to  conclude: 


a) 


b) 


Max.  G 


SNR 


(x  =  r^) 


|E°(0)|2 

a  -  V  cT1 


B 


(16) 


(17) 


Equations  (16)  and  (17)  give  the  results  of  a  first-order  perturbation. 

After  the  components  of  x^.  have  been  determined  from  (17) ,  they  can  be 
substituted  back  in  (9).  One  can  then  use  (D°  +  x^)  as  the  new  normalized 
element-position  column  matrix  and  perform  a  second-order  perturbation  to 
obtain  further  improvement  in  the  performance  index.  This  process  can  be 
repeated  until  it  becomes  evident  that  further  iteration  yields  a  negligible 
improvement.  For  the  many  cases  which  have  been  computed,  it  is  found  that 
convergence  toward  optimum  values  usually  takes  place  very  quickly;  seldom 
are  more  than  two  iterations  required. 
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The  spacing  perturbation  technique  outlined  above  yields  the  required 
element  positions  in  order  to  maximize  Gg^  for  a  given  set  of  excitation 
parameters.  Even  if  one  starts  with  a  uniformly  spaced  array,  the  element 
spacings  will  no  longer  be  uniform  after  the  perturbation.  Now  this  per¬ 
turbed,  nonuniformly  spaced  array  can  be  further  optimized  by  proper  ampli¬ 
fications  and  phase  shifts  following  the  array  elements.  With  this  in  mind, 
we  let  yo  "  I0»  yn  =  xn  cos  $n>  and  yN+n  =  Ir  sin  <{>n  (n=l,2, . . .  ,N) .  Further¬ 
more,  let  hQ4)  =  n/2,  hn(^)  =  cos  D^,  and  hN+n(^>  =  sin  D  ^(n-1,2 . N) . 

Equation  (5)  can  then  be  converted  to  the  following  form: 


where 


y"  h  h*  y 
r  _  o  0 

SNR  -,  =  -  * 

(18) 

y'  By 

?  =  ty„.  yr  y2,  •••.  y2N] 

(19) 

5o=[ho(0)>  hl(0)’  h2(0) . h2N(0)1 

(20) 

and  B  =  [b^.]  is  a  (2N  +  n)  *  (2N  +  n)  square  matrix  with 


2tt  it 


bij  =  hi  f  ^  f  hi^  hj(^)  w(9,<f>)  sin  0  d  0 


(21) 


0  0 


GSNR  in  is  exPressed  as  a  ratio  of  two  quadratic  forms.  Since  B 


2,3 


is  symmetric  and  positive  definite,  we  know  *  that  G  can  be  maximized  by 

oJNK 


choosing  the  column  matrix  y  =  y„: 

M 


yM  =  B_1  h 
M  o 


(22) 
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and  that 


Max.  G 


SNR 


(y  -  yM> 


h'  l"1  h 
o  o 


(23) 


Equation  (22)  completely  specifies  the  amplifications  and  phase  shifts 
required  to  obtain  the  maximum  possible  GgNR,  as  given  in  (23) ,  for  the 
perturbed  array.  If  the  performance  Index  of  interest  is  directivity 
<v  yield  th®  required  excitation  amplitudes  and  phases  to  maxi— 

mize  Go.  We  have  now  reached  a  second  submaximum,  which  may  possibly  be 
further  improved  by  holding  the  excitation  unchanged  and  again  perturbing 
the  spacings.  The  cycle  may  be  repeated  until  further  adjustments  are  no 
longer  worthwhile . 


It  is  to  be  emphasized  that  this  alternate  spacing  perturbation 
and  excitation  adjustment  procedure  of  seeking  a  maximum  performance  index 
can  be  applied  to  an  array  which  initially  has  an  arbitrary,  nonuniform 
spacing  and  an  arbitrary  distribution  of  excitation  amplitudes  and  phases. 
Optimization  using  the  above  technique  has  been  carried  out  for  both  broad¬ 
side  and  endfire  arrays  for  typical  spatial  distributions  of  noise  or  clutter 
power.  For  example,  with  a  cosine-square  over  a  pedestal  type  of  clutter 
power  distribution,  the  maximum  signal-to-noise  ratio  for  a  7-element  broad¬ 
side  uniform  array  is  16.0  and  that  for  an  optimized  array  is  increased  to 
181.9.  Even  more  imrtessive  improvement  is  obtained  for  a  7-element  endfire 
array,  in  which  case  the  signal-to-noise  ratio  is  increased  from  15.8  to  891  „ 4 
by  optimization.  Moreover,  the  improvement  does  not  result  in  the  objectionable 
features  (very  small  element  spacings  or  large  excitation  amplitudes  of  alter¬ 
nating  signs)  of  a  supergain  array. 
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III.  METHODS  FOR  OPTIMIZATION  OF  LARGE  ARRAYS 


Methods  for  the  optimization  of  either  the  directivity  or  the  signal- 
to-noise  ratio  of  an  N-element  antenna  array  requires  the  inversion  of  an 
N  x  N  square  matrix.  This  requirement  raises  questions  on  the  achievable 
accuracy  and  the  amount  of  computer  time  necessary  for  handling  arrays  with 
a  large  number  of  elements.  For  large  arrays,  it  would  be  extremely  desira¬ 
ble  if  a  technique  could  be  found  so  that  the  computational  difficulties 
would  be  reduced  in  the  use  of  the  optimization  method  based  on  a  considera¬ 
tion  of  the  ratio  of  two  Hermitian  forms.  Such  a  technique  has  been  found 
and  will  be  outlined  in  this  section. 

The  basis  of  this  technique  lies  in  a  transformation  which  transforms 
the  excitations'  in  the  elements  of  an  array  into  finite  summations  of  an 
exponential  series  with  complex  coefficients.  Let  the  normalized  array 
fact0*  wx  an  N-element  array  be  written  as 

N 

=  I  Ys  Xm  exp[j  (N  "  2m  +  *3  *  (24) 

m=l 

where  ip  has  been  defined  in  (3),  and  N  is  very  large.  We  set 

N 

Zm  55  ap  exp  [-j  (p-1)  (N  -  2m  +  1)  tt/N]  .  (25) 

P=1 

It  is  easy  to  verify  that  the  substitution  of  (25)  in  (24)  yields 

N 

BnW  =  i  Y.  apfpw  <26> 

P=1 
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(27) 


U*> 


sin  |  [if»  -  (p-1)  2ir/N] 
sin  |  [if;  -  (p-1)  2ir/N] 


Now  the  expression  in  (27)  has  the  following  property: 


where 


*(♦„  t)  =  N<5__t  » 

P  P  PP 


V  =  r  (p' ' 


(28) 


(29) 


and  6ppl  is  the  Kronecker  delta.  Combination  of  (27),  (28)  and  (29)  gives 
a  very  simple  result: 


E  (i|/  )  ■  a  , 
rivrpy  p  * 


(30) 


which  says  that  the  sampled  value  of  the  normalized  array  factor  in  the 
direction  specified  by  \p  is  numerically  equal  to  the  transformed  coefficient 
ap  in  (25).  Since  ap  decreases  rapidly  as  p  increases,  only  very  few  terms 
need  be  used  for  the  summation  in  (26) .  It  is  noted  that  the  inverse  rela¬ 
tion  for  (25)  is 


N 


-  |  £  Im  exp[j  (p-1)  (N  -  2m  +  1)  tt/N]  . 


(31) 


m=l 


For  a  properly  chosen  positive  integer  P  <<  N,  we  have  lap/ajj  «  1 


for  P  <_  p  <_  N-P,  and 


2P 


EnW  *  1 1]  Vp^  ’ 
P=1 


(32) 


where 


bp  =  ap  ,  gp(^)  =  fp(i/>)  for  1  1  P  1  P 


(33) 
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and 


bp  “  V2P+P*  8pW  *  £n-2p+pW  for  p  +  1  -  p  -  2P  •  (34) 

The  element  excitations  given  in  (25)  becomes 

P 

Im  ■  ^  bp  expi-j  (p-1)  (N  -  2m  +  1)  ir/N] 
p=l 

2P 

+  ^  bp  exp[-j  (N  -  2P  +  p  -  1)(N  -  2m  +  1)  it/N],  (35) 

p=P+l 

This  approximate  expression  for  I  can  then  be  used  in  the  formulation  for  the 

m 

array  performance  index  of  interest,  and  the  computation  labor  for  its  maximi¬ 
zation  is  greatly  reduced. 

Preliminary  investigation  has  indicated  that  8  values  of  give  a 
very  good  approximation  of  the  array  pattern  En(i|0  for  a  50-element  linear 
array.  This  means  that  the  tedious  and  inaccurate  inversion  of  a  50  *  50 
matrix  can  be  replaced  by  the  inversion  of  an  8  x  8  matrix  through  the  use  of 
transformation  (25) .  It  remains  to  obtain  more  numerical  data  for  arrays 
of  different  geometrical  configurations  and  for  different  scan  angles. 

Study  will  also  be  made  on  the  effect  of  the  approximation  on  meaningful 
constraints  such  as  main-beam  radiation  efficiency. 
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IV.  PATTERN  SYNTHESIS  FOR  CIRCULAR  ARRAYS  WITH  DIRECTIVE  ELEMENTS 


Circular  arrays  consisting  cf  radiating  sources  uniformly  disposed 
around  a  circular  ring  find  applications  in  situations  where  steerability 
and  pattern  invariance  are  of  importance.  Several  authors^-^  have  dealt 
with  pattern  synthesis  techniques  for  circular  arrays.  Most  published 
articles  pay  special  attention  to  patterns  with  equal  sidelobes  and  none 
appears  to  have  considered  nonisotropic  radiating  elements.  Unlike  a 
linear  array,  a  straightforward  multiplication  of  the  element  pattern 
function  and  the  array  factor  does  not  correctly  give  the  radiation  pattern 
of  a  circular  array.  This  causes  complication  in  pattern  analysis  and  dif¬ 
ficulty  in  pattern  synthesis .  It  should  be  pointed  out  that  the  problem 
of  pattern  synthesis  is  different  from  that  of  gain  optimization,  the  latter 
problem  having  been  solved  previously  for  circular  arrays  with  directive 
elements . ® 

o 

Recently  Mott  and  Dudgeon  proposed  a  method  for  synthesizing  the 
pattern  of  a  circular  array  with  isotropic  sources.  It  involves  a  prelimi¬ 
nary  integration,  leading  to  a  matrix  which  can  be  reduced  and  inverted  on 
a  digital  computer  with  less  labor.  However,  when  the  number  of  array 
elements  is  large,  error  accumulation  and  storage  problems  in  the  reduction 
and  inversion  of  a  large  matrix  would  again  become  serious.  In  this  section 
we  present  a  refined  method  which  does  not  require  a  matrix  inversion  at  all 
and  hence  is  particularly  useful  when  large  circular  arrays  with  many  elements 
are  to  be  designed.  Moreover,  the  present  method  is  applicable  to  cases  where 
the  array  elements  themselves  possess  a  directive  pattern. 
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Consider  a  circular  array  of  radius  a_  with  +  1  directive  elements 
arranged  uniformly  around  its  periphery.  If  g(ij>)  denotes  the  element  pattern, 
the  far-field  pattern  of  the  array  in  the  plane  of  the  circle  is 


M 

E(<l>)  =  L  Am8(4>  "  expU3a  cos  (*  “  fjf)]  »  (36) 

m=0 


where  3  =  2tt / X  and  Am  is  the  complex  excitation  in  the  mth  element.  In 
(36),  the  zeroth  element  is  assumed  to  lie  on  the  reference  axis  from  which 
the  angle  <j>  is  measured.  Let  Eq(<|>),  be  the  radiation  pattern  to  be  synthe¬ 
sized.  Since  Eq  (4>)  is  periodic  in  <f>,  it  can  be  expanded  in  a  Fourier-series 
form: 


where 


N 

Eo(40  -  V  Bn  exp(jn<j>)  , 
n=-N 


B 

n 


EQ(4>)  exp(-jn<J>)d(f>  . 


(37) 


(38) 


Equation  (27)  is  a  truncated  Fourier  series  of  2N  +  1  terms;  it  represents 
a  least-mean-square  approximation.  This  approximation  is  necessary  and 
reasonable  because  one  can  not  use  an  infinite  number  of  elements  in  a 
practical  design  and  because  the  avoidance  of  supergain  phenomenon  con¬ 
strains  the  minimum  allowable  element  spacing.  The  synthesis  problem  is 
then  the  determination  of  the  excitation  coefficients  {A  }  such  that  E(<M 

m  ' v  * 

in  (36)  equals  Eq (4>)  in  (37).  We  write 
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M  N 

L  Arns(*  -  m5  expIJBa  cos  <*  -  -  T  Bexp(jn*)  , 

m=0 


(39) 


n=-N 
for  0  <_  <f>  <_  2tt  . 


For  a  unique  solution  of  (39),  we  must  have  M  +  1  =  2N  +  1,  or 
M  ■  2N.  Integrating  both  sides  of  (39)  in  the  manner  of  (38),  we  obtain 

M 

Y  Am  °n  exPf“J2mnir/(M  +  D]  -  Bn  ,  (40) 

m=0 


where 


a 

n 


g(4>)  exp[j(8a  cos  <j>  -  n<j>)]d(j) 
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n-0,  +1,  +2,  ... 


+  N. 


(41) 


With  the  reasonable  assumption  that  the  element  pattern  g(<f>)  is  an  even 
function  of  <j>,  we  can  write 


K 

g(<J>)  = 


cos  kq>  . 


k=0 

Substitution  of  (42)  in  (41)  yields 


(42) 


a 

n 


Jn-k<6a>  +  J°+k 


^n+k 


(6a)]  , 


(43) 


n  =  0,  ±  i,  ±  2 . +  N  . 


Let  {Ck>  be  such  that  aR  *  0  for  all  n.  Then,  by  using  the  following 
identity, 
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N 

E 


exp[j2(m'  -  m)nir/(2N  +  1)]  =  (2N  + 


1)6  » 

mm 


(44) 


n»-K 

where  6^,  is  the  Kronecker  delta  (which  equals  zero  when  m  #  m'  and  equals 
unity  when  m  =  m’),  we  obtain,  from  (40), 


Am  "  liiVT  E  71  eXp[j  2mn,r/(2N  + 

r=-N  n 


(45) 


m  =  0,  1,  2,  . . ,  M(or  2N) . 

We  note  that  (45)  is  an  explicit  expression  for  finding  the  excitation 

coefficients,  (A  },  in  the  M  +  1  or  2N  +  1  elements  of  the  circular  array. 

m 

The  {Am>  so  determined  will  yield  a  radiation  pattern  which  approximates 
the  desired  Eq(4>)  in  the  least-mean-square  sense.  The  usefulness  of  (45) 
stems  not  only  from  its  applicability  to  the  synthesis  of  circular  arrays 
with  directive  elements,  but  also  from  the  fact  that  no  matrix  reduction  or 
inversion  is  required.  This  latter  feature  is  a  most  desirable  one  when 
large  circular  arrays  with  many  elements  are  to  be  synthesized. 
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V.  APPENDIX 


In  this  section  the  theorem  which  forms  the.  basis  of  the  spacing 
perturbation  technique  described  in  Section  II  is  stated  and  proved. 
Theorem.  If  a  quantity  P  can  be  expressed  in  terms  of  an  N  x  l  real  column 
vector  x  as 


P  =  A  -  2x'g  +  x’  C  x  ,  (46) 

where  A  is  a  constant,  $  is  another  N  *  1  real  column  vector,  x'  is  the 

I  a 

transpose  of  x,  and  C  is  an  N  *  N  positive-  definite,  symmetric,  square 
matrix,  then 

a)  min.  P  -  A  -  3'  C-1  3  ,  and  (47) 

(x  -  ^) 

b)  -  C'1  6.  (48) 

Proof:  If  C  is  positive  definite,  it  is  known  that10 

(S'C-1  B)(x'  Ci)>  (x’  I)2  (49) 

or 

x'  c  x  >~ —  (x?  S)2  ,  (50) 

S’  c  s 

where  the  equality  sign  applies  when 


Let  c  =  S' 


3  >  0, 


x  =  *M 
and  b  =  x 


*6.  We  have, 


from  (46)  and  (>0), 


But, 


P=A-2b+x'  Cx 
>  A  -  2b  +  — 


A  -  2b  +  —  =  A  -  c  +  —  (c-b)2 
c  c 

A  -  c  t 
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(51) 


(52) 


(53) 


Combining  (52)  and  (53),  we  obtain 


P  >_  A  -  c"1  3 

where  the  equality  sign  holds  with  (51);  hence  the  theorem  is  proved. 


(54) 


j 

| 

j 

\ 

j 

j 

.1 
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